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We study collective spin excitations of a magnetically ordered state in a multi- 
component system composed of both itinerant electrons and local moments. Here 
the induced spin-density- wave (SDW) ordering of itinerant electrons and the collinear 



i 

^ antiferromangetic (AF) ordering of local moments are locked together via a Hund's 

d; 

rule coupling. We show that the Goldstone theorem still holds at the RPA level with 
the gapless spin wave protected inside the small SDW gap of itinerant electrons, 

which, however, is fragile in the presence of ion-anisotropy. A gapped "out-of-phase" 

O ! 

spin mode extending over a much wider energy scale above the SDW gap is found 

to be more robust against the ion-anisotropy, which is mainly contributed by the 

local moment fluctuations. While the scattering between the Goldstone mode and 

itinerant electrons diminishes within the SDW gap, the "out-of-phase" mode will 



strongly interact with itinerant electrons and thus dominate the spin and charge 
dynamics in such an ordered phase. Possible relevance of such a model to the iron- 
pnictides will be also discussed. 

I. INTRODUCTION 

The recent discovery of high-T c superconductors in the iron pnictidesi has renewed a 
tremendous interest in the interplay between the magnetism and superconductivity. The 
similar issue has been vigorously investigated in the high-T c cuprates over two decades, where 
antiferromagnetism has been firmly related to the localized electrons of the Mott insulator 
in undoped cuprate compounds 2 -. By contrast, in the pnictides, the undoped parent material 
is not a simple Mott-Hubbard insulator but resembles a multi-band bad metal - i.e., the iron 



3<i-electrons are believed to be quite itinerant with their hybridized multi-orbitals forming 
multiple Fermi pockets at the Fermi level 3 - - -. 

It is natural for many to consider the spin-density-wave (SDW) order, observed in the 
undoped parent compounds by the neutron scattering measurements 8 , as originated from 
the same itinerant electrons via Fermi surface nesting- This picture seems more consistent 
in explaining the ARPES^ - — , transporti^, and optical properties 13 than a purely localized 
model, e.g., the J1-J2 modet^ - — , where a Mott-insulator transition has been implied. The 
latter is more reasonable in explaining the spin excitations in neutron scattering^ and 
high-temperature magnetic susceptibility^ - —, but obviously fails in understanding the bad 
metal behavior-^ 12 ^ 2 - and the presence of the small Fermi pockets^ - — i2&2£ in the parent 
compounds. 

Based on the overall experimental evidence, an alternative picture has been recently 
proposed^ 5 -, which assumes that some kind of orbital-selective Mott transition happens in 
the iron 3d orbitals of the pnictides such that both itinerant and Mott-localized electron 
coexist in the system. The minimal modepi based on this picture tries to reconcile the 
seemingly contradictory experimental facts and provides a natural understanding of the 
unified driving force behind the collinear antiferromangetic (AF) order and high-temperature 
superconductivity. Recently a microscopic realization of an orbital-selective Mott transition 
in the pnictides has been studied2&2£ based on the dynamic mean-field theory, which lends 
further support to this model. 

The key and unique feature for such a coexistent itinerant and localized electron system 
is that the two subsystems share the same characteristic momenta at Q s = (n,0) or (0, 7r). 
Namely, the hole and electron Fermi surface pockets around T and M points in the Brillouin 
zone (BZ) are approximately connected by Q s in the undoped case (i.e., close to the Fermi 
surface nesting), and at the same time, the local moments are strongly correlated at the 
AF wavevectors Q s . As a consequence, the local Hund's rule coupling between the itinerant 
electrons and local moments can be significantly enhanced around Q s , which is called the 
"resonant effect"—. At low temperature, such a "resonant effect" can serve as a predominant 
force in driving the magnetic or pairing instability at different dopings. Here the magnetic 
phase is predicted 25 to be an induced SDW order of the itinerant electrons locking with the 
collinear AF order of the local moments at the same Q s . Corresponding to such an AF 
ordering, a small SDW gap will open up in the excitation spectrum of itinerant electrons, 



although not necessarily pinned at the Fermi level as in an SDW state purely driven by 
Fermi surface nesting. At the mean-field level, the low-lying AF fluctuation of the local 
moments is also gapped at Q s due to the mutual locking of the magnetic orders in the 
two subsystems. Therefore, after the spontaneous magnetic symmetry breaking, the strong 
"resonant" scattering between the two degrees of freedom gets substantially reduced, which 
leads to a very coherent charge transport contributed by the ungapped part of the Fermi 
surfaces in consistency with the optical measurement^. 

However, it remains an important issue whether a gapless spin wave, i.e., a Goldstone 
mode, is still present in the magnetically ordered state of such a mult i- component system. 
To answer this question, one has to go beyond the mean-field theory to study the collective 
spin fluctuations, which is also important in order to self-consistently address the issue how 
the charge dynamics gets reshaped in the AF phase. In this paper, we shall address this 
issue with using a realistic five-band model^ to characterize the itinerant electrons near the 
Fermi pockets, and a J1-J2 type model to describe the Mott-localized electrons. Then we 
study the Hund's rule coupling between the itinerant and localized electrons at the RPA 
level in the magnetic ordered phase. We demonstrate that the Goldstone theorem indeed 
holds at the RPA level as a gapless spin wave emerges within the mean-field SDW gap 
of itinerant electrons. But it is fragile against the ion-anisotropy. We further find that 
the coupling between the Goldstone mode and the charge carriers diminishes in the long- 
wavelength around Q s as expected. On the other hand, distinct from a single component 
system, a gapped "out-of-phase" collective spin mode is also present with its high-energy 
part predominantly contributed by the local moments whose energy scale ~ Ji extends over 
a much wider regime than the SDW gap. Its low-energy part gets strongly renormalized 
by coupling to the itinerant electrons around Q s — it becomes gapped once the mean-field 
SDW order forms by itinerant electrons, which is not significantly modified at the RPA 
level. Furthermore, such an "out-of-phase" mode is not sensitive to a weak ion-anisotropy 
effect and is thus more robust than the Goldstone mode. The existence of the two branches 
of spin excitations in the AF state is a unique prediction of the present multi-component 
model. In particular, it is this "out-of-phase" spin mode that remains strongly interacting 
with itinerant electrons, at an energy higher than its gap, and therefore dominates the 
high-energy magnetic and transport properties in the magnetically ordered phase. 

The remainder of the paper is organized as follows. In Sec. II, we introduce the model 



and present the mean-field treatment in the magnetically ordered state. Then in Sec. Ill, 
we discuss the spin dynamics at the RPA level and demonstrate that the spin collective 
excitations are split into a gapless Goldstone mode which is upper-bounded and a gapped 
out-of-phase mode which extends over a much wider energy scale. In Sec. IV, we study the 
scattering between the collective spin modes and the itinerant electrons based on the single- 
particle self-energy of the itinerant electrons and the optical conductivity, which illustrate 
that the out-of-phase mode will play a dominant role beyond its energy gap. Finally, Sec. 
V is devoted to the discussion and conclusion. 

II. A MULTI-COMPONENT SYSTEM OF COUPLED ITINERANT AND 

LOCALIZED ELECTRONS 

A. Model 

We consider a multi-component system composed of coexistent multiband itinerant and 
Mott-localized electrons described by 

H = H it + H lo + H jH . (1) 

Here H lt is a tight-binding model of multiband itinerant electrons: 

-"it = — / , tij,mn c i m cT C jna, {*) 

i,j,m,n,cr 

where m and n are the orbital indices. The hopping integral £y, mn in H lt will be given based 
on a realistic five-band tight-binding model proposed 7 to describe the undoped iron-pnictide 
materials. The resulting band structure near the Fermi energy is shown in Fig. 1 by the 
solid (black) curves, with the Fermi surface shown in Fig. 2(a) in the undoped case (i.e., six 
electrons per site). Note that some slight modification with a global renormalization factor 
reducing the bandwidth has been phenomenologically made here in order to be consistent 
with ARPES^~— >2i. As shown in Fig. 2(a), the itinerant electrons form hole and electron 
pockets at the Fermi energy, which are located at the T point and the M point, respectively, 
separated by the momenta Q s in an extended BZ. 

The second term H\ Q in (pQ) describes the localized electrons in which a Mott gap is 
opened up via the so-called orbital-selective Mott transition 25 . Namely, the corresponding 



electrons only contribute to spin fluctuations, near the Fermi energy, by the local moments 
formed from the filled lower Hubbard band. Note that a microscopic realization of the 
orbital-selective Mott transition in such a system has been recently discussed based on the 
dynamic mean- field theory^ 1 ^. We shall simply use a J1-J2 model of S — 1 to depict the 
superexchange couplings between these local moments, i.e., 

H\ = J\ 2_^ Sj • Sj + J2 2_, Si • Sj (3) 

m am 

where (ij) and ((ij)) denote the nearest-neighbor and next-nearest-neighbor coupling, re- 
spectively. Here we assume Ji < 2J 2 such that the ground state of H\ Q itself may form 
a collinear AF ordering at the wavevector Q s . In principle, the same five-band electrons 
in iron-pnictides should contribute to both itinerant and local moment degrees of freedom. 
But for simplicity we treat H lt and H\ Q as if they govern independent degrees of freedom 
in the low-energy sector near the Fermi energy, so long as the Mott gap remains a large 
energy scale. We ignore the issues like how the Fermi surface shape gets affected by the 
orbital-selective Mott transition as well as how the Luttinger volume is correctly accounted 
for and mainly focus on the low-energy physics in the present work. 

The third term Hj H in (CQ) is the Hund's rule coupling between the spins of itinerant and 
the localized electrons: 

Hj H = — 2_^ J^^i " Sj m . (4) 

i,rn 

where s irn = c im aci m is the spin operator of the itinerant electrons in the m-th orbital and Sj 
denotes the localized moment at site i. J™ is a renormalized Hund's rule coupling constant. 
For simplicity, we shall assume a single J™ = Jq for different orbitals throughout the paper. 
Originally a simpler form of (CD) was proposed as a minimal model to describe the low- 
energy physics in the iron-pnictides^. The most important feature in such a model Hamil- 
tonian is that the peculiar momenta Q s , which on the one hand connects the two Fermi 
pockets of the itinerant electrons at V and M points, respectively, and on the other hand co- 
incide with the AF wavevectors of the local moments. It implies strongly enhanced dynamic 
coupling between the itinerant and local moment degrees of freedom, once the short-range 
AF correlations of the local moments set in around Q s even in a high-temperature normal 
state. With the decrease of temperature, such a "resonant" coupling will result in an AF 




FIG. 1: (Color online) The band structure near the Fermi energy for the five-band model described 
by Ha in © (black solid curves) in the undoped case. The reconstruction of the band structure 
of itinerant electrons in the presence of a collinear AF order is shown by red dotted curves. 

ordered phase with distinctive dynamic behaviors to be explored below, as compared to an 
ordinary SDW state of a pure itinerant electron system due to the Fermi-surface nesting 
mechanism or the collinear AF state of a pure J\-Ji mode. 

B. Mean-field treatment in the collinear AF ordered state 



In the following we first use a mean-field approximation to study the collinear AF ordered 
state in ([1]). The spin and charge dynamics at the RPA level will be investigated in the later 
sections. 

We start with the interaction term Hj H in (jlj) between the two sub-systems. By intro- 
ducing two order parameters of magnetization for the local-moments, M(i ), and itinerant- 
electrons, M( it ), respectively, as in \Sf) = M^ )e lCis ' r ' and J2 m (sf m ) = M^e lCis ' ri , one 
obtains the following linearization in Hj H , given by 



H Jh ->■ H i = -J 2 , 

i 
= #J(it) + Hj(i ). 



%^ s ' ri E^+%) eiQ> ' ri ^ 
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FIG. 2: (Color online) (a): The Fermi pockets of the itinerant electrons in the undoped case, (b): 
The reconstruction of the Fermi surface in the presence of a collinear AF order with a wavevector 
Q s = (vr,0) in the reduced BZ. 

Below the effect of such mean-field terms on the local moments and itinerant electrons will 
be explored in a self-consistent way. 



1. Local moment part 



In the magnetic order phase with (Sf) = M(i )e lC * 3 ' r % one may first use the conventional 
spin-wave approximation to treat H\ Q in (J3J) and then add .ffj(io) in (JHJ) to incorporate the 
effect of the coupling to the itinerant electrons. Here Q s is chosen to be (ir, 0). 



Introduce the Holstein-Primakoff (HP) transformation 



S 



LA 






2S - aldiCii, S iA = a[\/2S — a\a h S- A — S — a\a i: 



2S - b]b 3 b], Sj B = bJ2S-b% S Z JB = -S + b]bj, 



(6) 



where A and B sublattices are defined by the staggered factor e tCis ' ri = ±1. Under the 
approximation y 2S — a\ai ~ y 2S 1 — b\b{ ~ a/25* and using the boson operators in the 
momentum space 

1/2 



ai = ( jv) 2Z k akexp ( zk ' r 



H\ Q is transformed into 



H *> = S Yl t r k(4«k + b{b k ) + M^a^b^ + 4&L k )], 

^— 'k 

where ^ k means that the sum is within a reduced BZ with 

T k = 4 J 2 + 2Ji cos fcy, 
M k = 2J 1 cosk x + 4J 2 cosk x ■ cosky. 

Then ([9]) can be diagonalized by the Bogolubov transformation 

a k = u k a k + v k /3_\ k 
b\ = f k a fc + Mk/3l k , 



as follows 



where 



#io = 53 Wk(aiL«k + ^k), 



w k 



?'k 



w k 



^k + r k 



2w k 

-o; k + T k 

2w k 



^/r k -M^, 



sgn(M k 



(7) 
(8) 



(9) 



(10) 
(11) 



(12) 
(13) 

(14) 

(15) 
(16) 



Here o; k is gapless at k = Q s as the Goldstone mode of H\ Q in the AF ordered phase. 
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FIG. 3: (Color online) Dynamic spin susceptibility of local-moments at mean-field level, i.e. 
Xn ioq ( w ) = ~~ ^ m ^no)o ( < 1) < 1> W ) at q = Q s shows a gap, given by (fT9|) . which is opened up in 
the presence of a mean- field coupling to the SDW ordering of itinerant electrons, with the param- 
eters given in Sec. IIB. 

Now let us add #/(i ) in ([5]) arising from the mean- field decoupling of Hj H . It can be 
reexpressed in the spin-wave formalism by 



#/(io) = JoM(it) 5^ k ( a k a k + b l b k) + const. 



(17) 



This term will lead to a shift in I\ denned in (TTU1) : i.e., I\ — > I\ + JqMq^/S. As a result, 
the dispersion of the spin wave is modified by 



co k ^cu k = SJ(T k + J M (it) /S) - M, 



(18) 



In particular, u k is no longer gapless at Q s with an energy gap induced by M( it ) as 



u Qs = ^AM {it) SJ (2J 2 + J x ) + (M ( tt)Jo) • 
The order parameter M/y can be self-consistently calculated as 

M (lo) = S - i (^ (4a k + 6^ 



(19) 



(20) 
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where v^ is defined by (TT5l) with r k — > T k + JqM^/S and u^ defined in (TT8l) . 
Finally, the spin-spin correlations for the local-moment defined by 

xJo)(q,q,t) = -*(TS+(t)Sl q (0)), xJo)(q,q+Q s ,t) = -« (T5+(t)£ q _ Q .(0) 

can be obtained for this mean-field spin-wave state in the frequency space as 

1 



Xno)o( f i' q> w ) = M (i°) ( M k + ^k) 2 ( 



xJ o)0 (q>q + Qs> w ) = M (i°) 



+ i0+ -w - uj k + i0 + ' ' 

1 



oj — Uk + i0 + —to — Uk + i0+ 



(21) 



The spin gap in (TT9|) is clearly illustrated in the dynamic spin susceptibility x'n 



lo)0q 



W 



— 2Imxt^ (q, q, w) at q = Q s as shown in Fig. [3] (the parameters used are to be given 
below). It is pointed out that in the above calculation, we have further used the approxi- 
mation: S^/m —?• i v /2M(i )a] ( a/2M(i )6] ) in the original HP transformation ([6]) such that 
the spin commutation relations are satisfied at the mean-field level, i.e. (S x S) = ih(S 



2. Itinerant electron part 

Combining .ffj(it) in (jSJ) with the band kinetics energy term H lt in (}2"l) , the mean- field 
Hamiltonian of the itinerant-electrons reads 



Ha + #J(it) 



kmncr L 

JqM (Io) ct / t f \" 

„ mn I C k+Q s , mo- C k,ncr T c k,mcr C k+Q s ,n<7 I 



Z^ 



k(T 



^UffHffXiu,, 



where 



In the second line of 



Jmn V / / , v ij,mn^- 

Ti-Tj 

the 10 x 10 matrix i7 CT for the five bands is defined by 



(22) 



(23) 



(H) a 



with 



7 



F k -jla 
-lla F k+Qs 

J M M 



(24) 



(25) 
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Here I is the 5x5 identity matrix and F is the matrix defined by (F) m n = f mn . The column 
vector X kCT is given by 

^ko- — ( C k,laj C kj2o - • • • C k+ Q S) i cr C k+ Q S)2(T • • • ) • (26) 

By diagonalizing the 10 x 10 matrix (H) 

UlH a U a = D, (27) 

one gets 

#it + Hiat) — / „ -EkaC kQ , CT c ka(T , (28) 

where the band energy i? ka is equal to the a-th diagonal element of D, presented in Fig. 1 by 
the red dotted curves with the corresponding Fermi surface in the reduced BZ shown in Fig. 
2(b) in the undoped case. Here the order parameter M( it ) is self-consistently determined by 

in 

= JJ Tl Eka<EF E ^t (n + 5, a) £/kt (n, a) . (29) 

n 

Together with ([201), we find M (it) = 0.252 and M (lo) = 0.892 by choosing J = 20 meV at 
J\ = 0, and J 2 = 20 meV for S = 1. In the following we shall use these parameters to 
examine various spin and charge dynamics (we have also checked other small ratios of J1/J2 
and found the results remain qualitatively unchanged). 

Similar mean-field results have been previously obtained^ in a simpler model of (|TJ), where 
a spin gap similar to (fl9|) is also found in the spin-wave spectrum of local moments due to 
the Hund's rule coupling to the SDW order of the itinerant electrons. Such a gap will protect 
the collinear AF ordering jointly formed by both local moments and itinerant electrons below 
a transition temperature Tg^w- F° r example, the presence of this gap is reflected by a steep 
reduction of the uniform spin susceptibility below Tsdw, which is consistent 25 with the 
experimental measurement in the iron-pnictides^ 1 ^. However, according to the Goldstone 
theorem, a gapless mode is generally expected to exist in the AF ordered state. How such 
a Goldstone mode can be reconciled with the gapped local moment fluctuations discussed 
above will be studied at the RPA level in the next section. 

Finally, the spin-spin correlation functions of itinerant electrons are defined as 

Xjt7(q,q,*) = -*(TsJ(*)s: q (0)>, xfc(<l,<l + Qs,t) = -z(Ts+(t)5: q _ Q .(0)>, 
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FIG. 4: (Color online) Dynamic spin susceptibility of itinerant electrons at the mean-field level, 
i.e. X('it)oQ ( w ) — — 2ImX/j^Tn (Q»j Qs, 6 ^) with and without SDW order. A gap 2Asdw m the SDW 
ordered case is clearly shown. 



where s^ is the sum over all the five-orbital 



'? = £ 



'qm' 



(30) 



In the above mean-field state, after a straightforward but tedious calculation one obtains, 
for example, 

2 



X&)o W ^ 



E 



r(l) 



^ +q (a, /?) 



keft,k+qg,R, £ k , a >-E F , £ k+qi/3 <-E F 



W - (£ k ,a - ^k+q,/3) + «0+ 



+ 



E 



K 



(2) 



k,k+q±Q s 



:«,/?) 



kei?,k+q^fl, £ k , a >£ F , S k+q±Qsi/3 <E F 



W - (£k, a - £k+q±Q s ,/3) + ^ 



+ 



E 



Ki) 



^%,k 08, «: 



kefi,k-qefl, £ k , a >.E F , E k _ qi/5 <E F 



-w - (E Ka - -Ek-q,/?) + i0+ 



+ 



E 



K 



(2) 



k q±Q s ,k 



(£, a) 



keR.,k-q^iJ, £ k , a >£ F , S k _ q± q Si/3 <E F 



-a; — (Ek a — E) 



k-q±Q s 



z0- 



-,(31) 
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for the momentum q within the reduced BZ, i.e. q G R, and 

C (a, /?) = E Ki ( m > a ) ^t (m, /?) + U^ (m + 5, a) C/ qt (m + 5, /?)] 

m 

<q (a, = J2 i U H K «) ^t (™ + 5, /?) + t^ (m + 5, a) C/ qt (m, (3)} , (32) 

m 

with, say, 

5 =q= 5Z i.a^k+q.^k+q («>/?) + ^ ^.a^k+qiQ^t^k.k+qiQ, («> Z 3 ) • 

a^,ke_R,k+qe/J a/3,kei?,,k+q^i? 

(33) 
The other Green's function can be similarly obtained. Numerical calculation for the dy- 
namic spin susceptibility of the itinerant-electrons at the mean-field level, i.e. X(it)oQ ( u ) — 
— 2lmxtt) (Qs, Qs, w) is presented in Fig. 4 with and without SDW order, where the sup- 
pressing of the low-frequency spectrum by the SDW gap at 2A SDW is clearly shown. 

III. SPIN DYNAMICS 

In this section, we shall calculate the dynamic spin susceptibility at the RPA level be- 
yond the above mean-field approximation, by which both a Goldstone mode and a gapped 
collective spin mode in the AF ordered phase can be recovered. 

A. RPA treatment 

The Hund's rule interaction between the local moments and itinerant electrons in (j4j) is 
rewritten as 

Hj h = -JoY\ S q -s_ q (34) 

z — 'q 

where the spin operator § appears as a whole for the five-band. At the RPA level, for 
example, 

X+](q, q,t) = -i <TS+ (t) r q ) Q + l - <TS+ (t) H Jh fa) Hj H (t 2 ) s =q > + O (lf*J 

(T \ 2 »oo />0O 

y J / dt x J dtzlxpjo (q, q,t - h) X ^ )0 (q, q,h - t 2 ) X ^ )0 (q, q,* 2 ) 

+xJtH) (<i> q + Qs,* - *i) xjo) (q + Qs, q + Q s ,*i - h) xp )0 (q + Q s , q,^) 

+xj t y (<i> q + Q->* - *0 xjo) (q + Qs. q.*i - fe) xJ t 7o (q. q^)] + ..., (35) 
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FIG. 5: (Color online) Feynman diagrams of RPA for dynamic spin correlation of itinerant elec- 
trons. Here the single line bubbles denote the dynamic spin correlation function of (free) itinerant 
electrons, Xmq (Qj^) > the double line bubbles are the renormalized dynamic spin correlation func- 
tion of itinerant electrons, X(-*\ (q> w ) j the dotted lines are the dynamic spin correlation function 
of local moments, Xn^n ( < ^' u; )- The indices a, ft, label the matrix elements defined in (|36p . 



which is illustrated diagrammatically by Fig. |5j 

To make the formulation more compact, we define a 2 x 2 matrix x + ~ (q^) with the 
components 



xj7i) (q^) 



X 



X n 



x n 



X 



+- 



(q, q,^j 

(q,q + Q s ,u;) 

(q + Q s ,q + Q s ,w) 

(q+Q s ,q,cu). 



(36) 



With this definition, the Fourier transformation of ( |35l) reads 

X(it)(l,l) = X (it)0(l,l) + I y j 2^ X(it)0(l,i)X(lo)0(i,j)X(it)0(j,l) + 



*,i=l,2 



and more generally 



Jn 



X(it)(iJ) " = X('i t )0(i 1 i) + 1 2 ) ^ ^(it)0(i,i')^(lo)0(i',j')^(it)0(i',i) 4 



i'J'=l,2 



(37) 



(38) 
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FIG. 6: (Color online) Feynman diagrams of RPA for dynamic spin correlation of local moments. 
Here the double dotted lines are the renormalized dynamic spin correlation function of local- 
moments Xn~\ (q, ui). The other symbols are the same as in Fig. 5. 



As the Dyson equations in a compact matrix form, the above RPA result for the itinerant 

electrons can be reexpressed as 

-i 



xj it) (q, 



w 



J -(y) 2 Xct)o(q> w )xJ°)o(q> w ) 



X(it)o t q ' 



and similarly for the local moment part illustrated by Fig. |6l we have 

-l 



xJ o) (q,wj 



J, 



^2-+- 



J -(y) xj )o(q>w)xjt)o(q> 



u) 



xjo)o(q' w )' 



with the total dynamic spin susceptibility matrix given by 



xrpa (q> w ) = xj t ) (q> w ) + *(io) (q. w ) 



(39) 



(40) 



(41) 



B. Collective spin modes 



Now we focus on the total dynamic spin susceptibility defined by 

xrpa (q> w) = - 2Im XRP A (i,i) (q w ) ( 42 ) 

which can be numerically determined based on the RPA expressions (!39|) . (140]) and ( 14T|) 
given in the above subsection. 
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FIG. 7: (Color online) Total dynamic spin-susceptibility of coupled system in RPA level XrpaI 1 !' w ) 
fixing q = Q s . 

In Fig. [7J XrpaIQ) u ) a ^ a fixing q = Qs is shown as a function of u>. Here we find a sharp 
peak emerges at u — 0, in contrast to Figs. |3]andHl Namely a zero-mode (Goldstone mode) 
pole is indeed restored at q = Q s in the RPA spin-spin correlation function. Mathematically, 
it originates from the vanishing denominator in fl39|) and (j40l) . 

It is interesting to note that besides the Goldstone mode, there remains a high-energy 
mode as represented by the hump in Fig. [JJ It can be traced to the pole of Xn^o (*1> u ) ( c ^- 
Fig. [3]), only broadened through the scattering with the itinerant electrons. In other words, 
the gapped collective mode of the local moments identified in the mean-field state in the 
previous section is still present at the RPA level. Physically this gapped mode is an "out of 
phase" fluctuations of the local moments relative to the magnetization of itinerant electrons, 
while the Goldstone mode is the "in phase" fluctuations of the locked magnetizations from 
the local moments and itinerant electrons. 

To display the spin dynamics of the system in the whole BZ, the calculated Xrpa ( c Lj u> ) 
is presented in Fig. 8(a) with the x-axis representing the momentum q along the high- 
symmetry lines in the reduced BZ and the y-axis for the frequency u, while the brightness 
depicting the spectral weight Xrpa i.^ 00 )- From Fig. 8(a), it is clearly shown that the spin 
excitation spectrum in the AF ordered phase is separated into two branches, i.e. the lower 
Goldstone-mode branch and the upper "out-of-phase" mode branch. The dispersions of the 
two modes are illustrated in Fig. 8(b), which are defined by (q, oS) at the largest (brightest) 
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FIG. 8: (Color online) (a) The brightness represents the spectral weight of the dynamic spin 
susceptibility Xrpa m ^ ne 1 an< ^ w space, (b) The dispersions of two branches of the collective spin 
mode as read from (a). The dispersions of the spin- wave of the pure J\ — J<i model (dotted) and 
that with a gap opening at the mean-field level (dashed) are also shown for comparison. 



Xrpa ( f l5 a; ) m Fig- 8(a). The spin wave dispersion of the pure J\ — J 2 model, i.e., ( 1TB]) and 
that of the gapped one, (ITS]) , due to coupling to the SDW order of the itinerant electrons 
are also shown for comparison. 

The lower branch Goldstone mode in Fig. 8(b) is well-defined within a small region 
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which centers at the AF wavevector Q s , upper-bounded by the SDW gap, i.e. 2A SDW 
of the itinerant electrons, beyond which it decays quickly due to the scattering with the 
particle-hole continuum of itinerant electrons (cf. Fig. Hj). Inside the SDW order gap, such 
a Goldstone mode is protected and is decoupled from the itinerant electrons due to the 
reconstruction of the Fermi surface shown in Figs. Q] and EJ 

On the other hand, the high energy "out-of-phase" mode is mainly contributed by the 
local moments and present throughout the BZ. It gets slight renormalization and broadening 
due to the scattering from the itinerant electrons at the RPA level, but more or less follows 
the dispersion obtained in the mean- field treatment (dotted), with a gap r\ 

ri~u Qa . (43) 

It is noted that the Goldstone mode will be fragile against the ion-anisotropy. By adding 
an ion-anisotropy term, 

tf(ioB) = ~J Z E (^) 2 ( 44 ) 

i 

the spin- rotational symmetry will be broken in H\ Q . By using the spin- wave expansion 

#(ion) = 25 J z ^^ ( a k a k + ^k^k) + const. (45) 

one finds the spin-excitation spectrum (fl6|) acquires an ion-anisotropy gap 



w k = ^(r k + 2J z r-M2. (46) 

At the RPA level, the Goldstone spin wave will also become gapped and its spectrum weight 
is reduced with reference to that without ion-anisotropy, as clearly shown in Fig. [9j When 
one further increases J z , the Goldstone mode disappears, but the gapped "out-of-phase" 
collective mode remains robust and insensitive to the ion-anisotropy. 

Therefore, in the present coupled local moment and itinerant electron system, the Gold- 
stone theorem still holds for a spontaneously symmetry breaking state with AF ordering 
at the RPA level. However, the Goldstone mode is very sensitive to the presence of ion- 
anisotropy and does not play an important role for charge dynamics as it is decoupled from 
the itinerant electrons. On the other hand, the gapped "out-of-phase" collective mode is 
more prominent which extends over the whole BZ with a spin-wave bandwidth ~ 4J2 and 
can be easily probed by neutron scattering experiments. Such two-branch collective spin 
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FIG. 9: (Color online) The dynamic spin susceptibility similar to Fig. El but with an additional 
ion-anisotropy J z = 1 meV. 

excitations are unique prediction for the AF ordered phase. In the next section, the charge 
response will be further examined based on the scattering of such collective spin modes with 
the itinerant electrons in SDW ordering. 

IV. CHARGE DYNAMICS 

In this section, we study the charge dynamics in the coupled local-moment and itinerant- 
electron system in the AF ordered phase. Here the itinerant electrons in the Fermi surface 
region will be highly coherent even in the presence of the gapless Goldstone mode. 
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A. Self-energy of itinerant electrons 

By going beyond the mean-field linearization in (JS)) , we shall consider the scattering based 



on 



H Jh -». H'j = --^= Yl c L c ^t [ S U V & («. P) + 5 J-k+QXq («> P) 

V kei?,q£-R,a^ 



+ h.c. 



- -r4= E c U^t<^ + ^- (47) 

ZViV keR,qe-R,a/3 

/I Q\ 

where the V^ ' (ct,/3) are defined by (E2D- Here we retain only the scattering of itinerant 
electrons with the transverse fluctuations of local moments as the longitudinal fluctuations 
in S z are gapped with (S z ) ^ 0. 

The single-particle Green's Function can be evaluated perturbatively by 

G a0i (k, u ) = G\ (k, w) + E ^(k, u;)E 74 (k, w)G^(k, w) + O (iff) , (48) 
where the self-energy 

E 74 (k, W ) = ^ E y _ — g^ (q, u - n) ^ kq7 ^ 6 ( fi ) > ( 49 ) 



qGfl,^ " " x 



with 



*&„*«(«) = ^(^(q-k^)!^! +xJ o)(li2) (q-k ! fi)V«0 2 )* + 

x^ (2>1) (q - k, n) f^ 1 )* + x^ (2 , 2) (q - k, n) |y( 2 ) | 2 . (50) 



Here the matrix Xn^\ is the dynamic spin correlation function for the local-moment in the 
RPA level defined in (|40|) . Fig. 10 and Fig. 11 show the corresponding Feynman diagrams 
of the self-energy correction for itinerant electrons. 

The equation (jlSjl can be further expressed in a 10 x 10 matrix formalism by 

= G kiV + Gk^Sk^Gk^ (51) 

which gives rise to G^^ = ( (G k)CJ ) — (£k,w) ) • Here we have omitted the spin index as 
they are in fact spin- independent. 
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FIG. 10: (Color online) Feynman diagrams of self-energy correction for itinerant-electrons. Here 
the single black lines denote the proporgators of (free) itinerant-electrons, G°r(k, w), the double 
black lines are the renormalized proporgators of itinerant-electrons, G a g(k.,cj), the dotted wavy 
lines are the modified dynamic spin correlation of local-moments in the RPA level xTmir ac (^) 
defined by ([50]) . The indices a, /?, 7, p, ?, label different bands. 



By noting that the zero-th order single-particle Green's function is diagonal, through the 
Lehmann representation we have for ui > 

72 /-+oo jrj 

- ImS 7e (k, u) = ^J2 T P ° 6 (q ' U ~ Q) D S ( k " <*' fi ) ^ ^ + r^ (fi - w)] , 

(52) 

(53) 
(54) 



where 



and 



9°(q,a;) = -2ImG° e (q,u;) 



D+;(k-q,n) = -2Imx 



(i?)kq76»ee 



(n) 



are the spectral functions of single-particle and the modified dynamic spin susceptibility of 
the local-moment, respectively. 

At zero temperature, the imaginary self-energy (152J1 can be further evaluated as (for 
u > 0) 



— ImE Q!j a(k, u 



8N 



Yl Kep^-^^-E^). 



(55) 



qeB,,e,o<E liS <u 



2U(*,<») = 



%(R)kqr0& ^ ' 



G^faa-O) 
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FIG. 11: (Color online) Self-energy correction for itinerant-electrons 
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FIG. 12: (Color online) The imaginary part of the quasiparticle self-energy -ImS QQ ,(k,o;) at two 
typical k-points marked in Fig. 2(a) is shown. The red-circled curve is for incorporating the 
scattering with the full collective mode at the RPA level and the black-squared one is for the 
scattering with the bare local moment fluctuations governed by H\ . 



In Fig. 12, the -ImE aQ ,(k, ui) as a function of a; is shown at two typical momenta, (ki, a±) 
and (k 2 ,«2), respectively, which are marked at Fermi pockets in Fig. 2(a). It shows that 
the life time of the quasiparticle excitations at these points of Fermi pockets actually gets 
substantially enhanced at low u in the collinear AF ordered state (solid circles), as compared 
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to an artificial case (solid squares) without an induced SDW order appearing in itinerant 
electrons such that there is no gap rj opened up in the "out-of-phase" mode. It implies 
that although in the normal state the itinerant electrons may be strongly scattered by the 
low-lying local moment fluctuations, the sharp coherence of quasiparticles will emerge in the 
AF state, where the gapless Goldstone mode is essentially decoupled from the particle-hole 
continuum and the "out-of-phase" mode is gapped. 

B. Optical conductivity 

Finally, let us examine the overall structure of the optical conductivity a (to) at the mean- 
field and RPA levels. It is related to the current- current correlation function Gj (u) through 
the relation 

ff(w) = --ImGj(w), (56) 

CO 

where Gj (to) is the Fourier transformation of the current-current correlator 

Gj(t) = -i(TJ(t)J(0)). (57) 

Here J denotes the q = current operator J along, say, the x axis for the five-band model 
(ED defined by 

_ dH it {A) i 

= 2_^ Vk/mn(kjC kmcr Ckno- 
kcrmn 

= ^2 c laa c W* Vkfmn(k)U£ a (m,a)U ka (n,/3) + ^ k f mn (k+Q s )U£ a (m + 5,a)U ka (n + 5,l3) 

k£R,crmnaf3 

= Yl v *™pcL<T c w0- ( 58 ) 

keR,aa(3 

Omitting the vertex-correction, we find 
Gj («) = -^ E / ^V*»*V*«Jt/G(fa* ( k > fi ) <V„ (k, n + u). (59) 

with the Feynman diagrams shown in Fig. 13. Denoting 

p a/3 (k, w) = -2ImG a p a (k, u) sgn (u) , (60) 



V, 



k after 



G ' a jk,m 



G Pa , a (i,n+J) 
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FIG. 13: (Color online) Current-current correlation function for itinerant-electrons. Here the 
double black lines are the renormalized proporgators of itinerant-electrons. 



at T = we obtain 



1 f° dQ 

ImGj = — J2 J ^^m/^k^yp^'jk,!)) 



j(k,n 



to 



(61) 



kg_R, a,a/3,a 



such that 



a(u) 



1 



Nu 



£ 



dVL 

2% 



VtorapVtoa'ffPifa (K ^) Pfi a ' (K ^ 



U> 



(62) 



k£R,a,a0,a / 

The calculated optical conductivity a [uj) is shown in Fig. 14. The black-squared curve 
shows the result for the bare five-band itinerant electrons and the red circles represent 
the result of the SDW reconstructed bands. Note that the multi-peak structure is mainly 
due to the multi-band effect with the lower ones change significantly in the SDW state. 
Furthermore, by incorporating the scattering with the collective spin modes at the RPA 
level, as shown by the blue triangles in Fig. 14, no significant change has been found in the 
optical conductivity with q = 0. It clearly illustrated that the itinerant electrons remain very 
coherent in the AF ordered phase, where the low-lying Goldstone mode does not strongly 
scatter the quasiparticles as expected. 
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FIG. 14: (Color online) The optical conductivity a (w) in the AF ordered state (blue solid triangle) 
which shows that the quasiparticles remain quite coherent as compared to the case without incor- 
porating the scattering with the collective spin modes (red solid circles). For comparison, (t(uj) for 
the pure five-band model without the SDW reconstruction is shown (black square). 

V. DISCUSSION AND CONCLUSION 



In this work, we have studied the collective spin excitations in the AF ordered phase of 
a multi-component system composed of coexistent itinerant and localized electrons. The 
main prediction is that a usual spin mode is split into two branches in such a multi-band 
system with orbital-selective Mott transition. The lower branch is a gapless Goldstone mode 
which is recovered at the RPA level and is quickly damped above 2A SDW by coupling to 
the particle-hole continuum of itinerant electrons, similar to the case in which the SDW 
order is due to the pure Fermi-surface nesting effect for itinerant electrons. However, an 
upper branch remnant spin wave reemerges above the SDW gap over a much wider energy 
~ J2 which is dominantly contributed by the local moment fluctuations. Here the lower 
and upper branches can be regarded as in-phase and out-of-phase combinations of the spin 
fluctuations from the itinerant and Mott-localized electrons, which are clearly distinguished 
from the single mode in a conventional AF state either due to the pure Fermi-surface nesting 
effect for itinerant electrons or AF superexchanges of local moments. 
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Experimentally the high-energy spin-wave excitation has been clearly observed by the 
neutron scattering experiments over an energy scale ~ J2 and presumably survives in 
the high-temperature regime above the ordered phase. However, a small gap (~ 6 — 10 
meV) has been generally found in SrFe2As2 and BaFe2As2, and interpreted as due to ion- 
anisotropy^ 1 ^. As shown in this work, a small ion-anisotropy can indeed easily destroy the 
lower branch Goldstone mode, while the upper branch is more robust. It remains to be 
seen if the lower branch spin mode can be unambiguously identified for a sample with less 
ion-anisotropy. 

Another distinct property of the present system is that the itinerant electrons become 
very coherent in the AF ordered phase, leading to a good metallic behavior after the AF 
transition. This is in contrast to the presumably strong scattering between the itinerant 
electrons and local moments in the normal state, where due to the very fact that the mo- 
mentum displacement of the hole-electron Fermi pockets of the itinerant electrons matches 
with the AF wavevector of the local moments, there exists a strongly enhanced interaction, 
i.e., the "resonant effect" around Q s , between the two subsystems. It provides the strong 
scattering source responsible for a drastic change in the charge response once the system 
enters the AF long-range ordered state at low-temperature. By forming a joint magnetic 
ordering, the two subsystem effectively get "decoupled" as the AF fluctuations of the local 
moments gain a gap i] as the out-of-phase collective mode. On the other hand, the gapless 
Goldstone mode is effectively decoupled from the itinerant electrons in the long-wavelength, 
thanks to the Fermi surface reconstruction by the SDW order. 

Therefore, the collective fluctuations of the local moments can serve as the main driving 
force for both the AF ordering as well as the superconducting pairing in the system via the 
"resonant effect" on itinerant electrons, as first pointed out in Ref.— . The resulting magnetic 
and charge properties in the AF ordered state, in particular the two branch collective spin 
modes predicted in the present work, can be further tested by experiment in order to establish 
the relevance of the model with the iron pnictides. 
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